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1.  Euler  Integration  in  o-minimal  setting 

1.1.  Introduction.  Integration  with  respect  to  Euler  characteris¬ 
tic  is  a  homomorphism  fy  -dx  ■  CF{X)  TL  from  the  ring  of  con- 
structible  functions  CF{X)  ("tame"  integer-valued  functions  on  a 
topological  space  A")  to  the  integers  Z.  It  is  a  topological  integration 
theory  which  uses  as  a  measure  the  venerable  Euler  characteristic 
X-  Euler  characteristic,  suitably  defined,  satisfies  the  fundamental 
property  of  a  measure: 

(1.1)  x(^  U  5)  =  x{A)  +  xiB)  -  xiA  n  B).  ^ 

for  A  and  B  "tame"  subsets  of  X.  We  extend  the  theory  to  R-valued 
integrands  and  demonstrate  its  utility  in  managing  incomplete  data 
in,  e.g.,  sensor  networks. 

1.2.  Constructible  integrands.  Because  the  Euler  characteristic 
is  only  finitely  additive,  one  must  continually  invoke  the  word  "tame" 
to  ensure  that  x  is  well-defined.  One  means  by  which  to  do  so  it  via 
an  O-MINIMAL  STRUCTURE  [39],  a  sequence  O  =  (C>„)  of  Boolean  al¬ 
gebras  of  subsets  of  R"  satisfying  a  small  list  of  axioms:  closure  un¬ 
der  products,  closure  under  projections,  and  finite  decompositions 
in  Oi-  Elements  of  O  are  called  DEFINABLE  sets  and  these  are  "tame" 
for  purposes  of  integration  theory.  Examples  of  o-minimal  struc¬ 
tures  include  (1)  piecewise-linear  sets;*  (2)  semi-algebraic  sets;  and 
(3)  globally  subanalytic  sets. 

Definable  functions  between  spaces  are  those  whose  graphs  are 
in  O.  For  A'  and  Y  definable  spaces,  let  Def(A'.  I')  denote  the  class 
of  compactly  supported  definable  functions  fi  :  X  V,  and  fix 
as  a  convention  Def(A^)  =  Def(A',  R).  Let  CF(A')  =  Def(.V,Z)  C 
Def  (A',  R)  denote  the  ring  of  CONSTRUCTIBLE  FUNCTIONS:  compactly 
supported  Z-valued  functions  all  of  whose  level  sets  are  definable. 
Note  that  in  general,  definable  functions  (even  definable  'homeo- 
morphisms')  are  not  necessarily  continuous. 

We  briefly  recall  the  theory  of  Euler  integration,  established  as 
an  integration  theory  in  the  constructible  setting  in  [28,  36,  35,  38] 
and  anticipated  by  a  combinatorial  version  in  [6,  21,  22,  34].  Fix  an 
o-minimal  structure  O  on  a  space  A'.  The  geometric  Euler  character¬ 
istic  is  the  function  y  :(!?—»  Z  which  takes  a  definable  set  .4  C  O  to 
X{A)  =  X^,(— 1)' dim R),  where  is  the  Borel-Moore  ho¬ 
mology  (equivalently,  singular  compactly  supported  cohomology) 

*Some  authors  require  an  o-minimal  structure  to  contain  algebraic  cur\-es, 
eliminating  this  particular  example. 
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of  A.  This  also  has  a  combinatorial  definition;  if  A  is  definably  home- 
omorphic  to  a  finite  disjoint  union  of  (open)  simplices  ]J^  cTj,  then 
x(j4)  =  — Algebraic  topology  asserts  that  \  is  indepen¬ 

dent  of  the  decomposition  into  simplices.  The  Mayer-Vietoris  prin¬ 
ciple  asserts  that  x  is  a  measure  (or  'valuation')  on  O,  as  expressed 
in  Eqn.  [1.1]. 


The  Euler  integral  is  the  pushforward  of  the  trivial  map  A' 
{pt}  to  fy  dx  :  CF{X)  — >  CF{{pt})  =  Z  satisfying  Jy  1^  d\  =  \(A) 


for  1^  the  characteristic  function  over  a  definable  set  A.  From  the 
definitions  and  a  telescoping  sum  one  easily  obtains: 


oo 


(1.2)  \  hd\  =  sx{h  =  s}  =  '^xih  >  <  -■‘’I- 

\  c—  —  -v-  c — n 


Because  the  Euler  integral  is  a  pushforward,  any  definable  map  F 
X  induces  F.  :  CF{X)  CF{Y)  satisfying  h  dx  =  fy  FJi  d\ . 
Explicitly, 


I 


FMy)  = 


(1.3) 


hd\. 


as  a  manifestation  of  the  Fubini  Theorem. 

The  Euler  integral  has  been  found  to  be  an  elegant  and  useful  tool 
for  explaining  properties  of  algebraic  curves  [8]  and  stratified  Morse 
theory  [37, 9],  for  reconstructing  objects  in  integral  geometry  [35],  for 
target  counting  in  sensor  networks  [1],  and  as  an  intuitive  basis  for 
the  more  general  theory  of  motivic  integration  [12, 13]. 

1.3.  Real-valued  integrands.  We  extend  the  definition  of  Euler 
integration  to  R-valued  integrands  in  Def(A')  via  step-function  ap¬ 
proximations. 

1.3.1.  A  Rieniann-sum  definition. 

Definition  1.  Given  h  e  Def(A"),  define; 


(1.4) 


(1.5) 


We  establish  that  these  limits  exist  and  are  well-defined,  though 
not  equal. 


Lemma  2.  Given  an  affine  function  h  G  Def(a)  on  an  open  k-siniplex 

o, 

(1.6)  f  /t  [c/xj  =  (~1)^  iiif /i ;  f  /t  ffi\]  =  (-!)*■’ Hup/i. 

J  (j  ^  J  a 

This  integration  theory  is  robust  to  changes  in  coordinates. 

Lemma  3.  Integration  on  Def(A^)  with  respect  to  [dxj 
invariant  under  the  right  action  of  definable  bijections  of  X. 

Lemma  4.  The  limits  in  Definition  1  are  well-defitted. 

Integrals  with  respect  to  [dxj  ^nd  [dxl  ^re  related  to  total  varia¬ 
tion  (in  the  case  of  compactly  supported  continuous  functions). 

Corollary  5.  If  M  is  a  1-dimensional  manifold  and  h  G  Def(3/) 
is  continuous,  then 

(1.7)  [  h  [dxj  ^  -  [  h  [dvl  =  ^totvar(/i). 

Jm  Jm  ^ 

This  result  generalizes  greatly  via  Morse  theory;  see  Corollary  11. 
One  notes  that  [dx\  and  \dx]  give  integrals  which  are  conjugate  in 
the  following  sense. 


Lemma  6. 

(1.8)  J  h\d\]  =  -  J -h[d\\. 

The  temptation  to  cancel  the  negatives  must  be  resisted:  see  Lemma 
12  below. 

1.3.2.  Computation.  Definition  1  has  a  Riemann-sum  flavor  which 
extends  to  certain  computational  formulae.  The  following  is  a  defin¬ 
able  analogue  of  Eqn.  [1.2]. 


Proposition  7.  For  h  e  Def (A'), 


(1.9)  [  h  [dxj  =  [  x{h  >  s}  -  \  {h  < 

Jx  Js-{) 

(1.10)  /  /i[dxl  =  r  x{h  >  s}  -  x{h  < 

Jx  Js^O 

It  is  not  true  that  h  [dxj  =  =  s}ds: 

generalization  of  Eqn.  [1.2]  is  the  following: 


—  .s  }  ds 
— ,s}  ds. 

the  proper  Lebesgue 
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Propositions.  For  he  Def(A'), 

(1.11) 

/  h  [(i\J  =  liin  -  [  sx{ 
Jx  '—0+  (  7r 

(1.12) 

[  h  [d\']  =  liin  -  [  s  \  { 

Jx  f  Jr 

<  h 

<  h 


<  6-  +  t}  ds 

<  ,s  +  f}  ds. 


133.  Morse  theory.  One  important  indication  that  the  definition 
of  f  [d\J  is  correct  for  our  purposes  is  the  natural  relation  to  Morse 
theory:  the  integrals  with  respect  to  [dx\  and  [d\]  are  Morse  index 
weighted  sums  of  critical  values  of  the  integrand.  This  is  a  localiza¬ 
tion  result,  reducing  from  an  integral  over  all  of  A’'  to  an  integral  over 
an  often  discrete  set  of  critical  points. 

Recall  that  a  function  h  :  M  R  on  a  smooth  manifold  M 
is  Morse  if  all  critical  points  of  h  are  nondegenerate,  in  the  sense 
of  having  a  nondegenerate  Hessian  matrix  of  second  partial  deriva¬ 
tives.  Denote  by  C{h)  the  set  of  critical  points  of  h.  For  each  p  G  C(//), 
the  Morse  index  of  p,  fi{p),  is  defined  as  the  number  of  negative 
eigenvalues  of  the  Hessian  at  p,  or,  equivalently,  the  dimension  of 
the  unstable  manifold  W^{p)  of  the  vector  field  —V/?  at  p. 

Stratified  Morse  theory  [20]  is  a  powerful  generalization  to  trian- 
gulable  spaces,  including  definable  sets  with  respect  to  an  o-minimal 
structure  [9,  37].  We  may  interpret  [dx]  and  \d\]  in  terms  of  the 
weighted  stratified  Morse  index  of  the  graph  of  the  integrand. 


Definition  9.  For  A"  C  R^*^  definable  and  h  G  Def(A'),  cdefine  the 
co-index  of  h,  X*h  to  be  the  stratified  Morse  index  of  the  graph  of  h, 
F/,  C  A"  X  R,  with  respect  to  the  projection  tt  :  A"  x  R  — >  R: 

(1.13)  {X*h){x)  -  Hill  Y  n  {h  <  h{x)  +  f'}")  , 


where  B^{x)  is  the  closed  ball  of  radius  t  about  x  G  A".  The  index  I, 
is  the  stratified  Morse  index  with  respect  to  height  function  —  tt:  i.e., 
XJi  =  or 

(1.14)  {X^h)ix)  -  liin  V  (b,{x)  H  {h  >  h{x)  -  t'})  . 


Note  that  :  Def(A^)  — ^  CF{X),  and  the  restriction  of  these 

operators  to  CF{X)  is  the  identity  (every  point  of  a  constructible 
function  is  a  critical  point).  The  two  types  of  integration  on  Def(A') 
correspond  to  the  Morse  indices  of  the  graph  with  respect  to  the  two 
orientations  of  the  graph  axis  —  the  projections  tt  and  —  tt. 
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Theorem  10.  Tor  any  continuous  h  G  Def(A"), 

(1.15)  J  h  [f/xj  =  hr  h  dx  ;  j  h\  dx  ]  =  ^  hi.  hd\. 

Corollary  11.  If  /i  is  a  Morse  function  on  a  closed  n-manifold 


M,  then; 

(1.16) 

/  h[dx\=  Y. 

p€C(h) 

(1.17) 

pec{h) 

From  this,  one  sees  clearly  that  the  relationship  between  [r/x  J  and 
[dy]  is  regulated  by  Poincare  duality.  For  example,  on  continuous 
definable  integrands  over  an  n-dimensional  manifold  M , 

(1.18)  /  A  (r/xl  =  (-!)"  /  h[dx\. 

Jm  jm 

The  generalization  from  continuous  to  general  definable  integrands 
is  simple,  but  requires  weighting  l*}i  by  h  directly.  To  compute 
Jy  h[dx\f  one  integrates  the  weighted  co-index 

(1.19)  lim  h{x  +  f')\  n  {h  <  h{x)  +  t'}) 

with  respect  to  r/y. 

Corollary  11  can  also  be  proved  directly  using  classical  handle- 
addition  techniques  or  in  terms  of  the  Morse  complex,  using  the 
fact  that  the  restriction  of  the  integrand  to  each  unstable  manifold 
of  each  critical  point  is  unimodal  with  a  unique  maximum  at  the 
critical  point.  It  is  also  possible  to  express  the  stratified  Morse  index 

—  and  thus  the  integral  here  considered  —  in  terms  of  integration 
against  a  characteristic  cycle,  cf.  [20,  37]. 

One  final  means  of  illustrating  Corollary  11  is  to  use  a  deforma¬ 
tion  argument.  Let  h  be  smooth  on  and  (ht  be  the  flow  of  —V/?. 
Then  the  integral  is  invariant  under  the  action  of  (pt  on  h;  yet  the  lim¬ 
iting  function  —  lim^-^oc  h  o  is  constant  on  stable  manifolds  of 

-  V/i  with  values  equal  to  the  critical  values  of  li.  We  have  not  shown 
that  the  limiting  function  is  constructible  (this  depends  on  the  exis¬ 
tence  of  definable  invariant  manifolds  —  we  are  unaware  of  relevant 
results  in  the  literature)  and  thus  do  not  rely  on  this  method  for  proof 
but  rather  illumination. 
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1.4.  The  integral  operator.  We  consider  properties  of  the  inte¬ 
gral  opera tor(s)  on  Def(A'). 

1.4.1.  Linearity.  One  is  tempted  to  apply  all  the  standard  con¬ 
structions  of  sheaf  theory  (as  in  [36,  35])  to  f  y  :  Def(A')  R.  How¬ 
ever,  our  formulation  of  the  integral  on  Def(A)  has  a  glaring  disad¬ 
vantage. 

Lemma  12.  fy  :  Def(A)  ^  R  (via  [riyj  or  [dx])  is  not  a  homomor¬ 
phism  for  dim  X  >  0. 

This  loss  of  functoriality  can  be  seen  as  due  to  the  fact  that  [/  +  r/J 
agrees  with  [/J  -t-  [gj  only  up  to  a  set  of  Lebesgue  measure  zero,  not 
y-measure  zero.  The  nonlinear  nature  of  the  integral  is  also  clear 
from  Eqn.  [1.15],  as  Morse  data  is  non-additive. 

1.4.2.  The  Fiibini  Theorem.  In  one  sense,  the  change  of  variables 
formula  trivializes  (Lemma  3).  The  more  general  change  of  variables 
formula  encapsulated  in  the  Fubini  theorem  does  not,  however,  hold 
for  non-constructible  integrands. 

Corollary  13.  The  Fubini  theorem  fails  on  Def(A')  in  general. 

Fubini  holds  when  the  map  respects  fibers. 

Theorem  14.  For  h  e  Def(A'),  let  F  :  X  Y  be  definable  and  h- 
preserving  (h  is  constant  on  fibers  of  F).  Then  fy  F*/;  [r/\J  =  fy  /)  ld  \  ], 
and  Jy  F,h\d\]  =  /i[dx]. 

Corollary  15.  For  h  €  Def(A'),  jy  h  =  hji.  In  other  words, 


(1.20) 


and  likewise  for  [dy]. 

1.4.3.  Continuity.  Though  the  integral  operator  is  not  linear  on 
Def(A),  it  does  retain  some  nice  properties.  All  properties  below 
stated  for  f  [d\J  hold  for  f  [d\']  via  duality. 

Lemma  16.  The  integral  f  [dyj  :  Def(A')  — >  R  is  positively  homoge¬ 
neous. 

Integration  is  not  continuous  on  Def(A")  with  respect  to  the  C‘' 
topology.  An  arbitrarily  large  change  in  j  /i[dxj  may  be  effected  by 
small  changes  to  h  on  a  (large)  finite  point  set.  In  some  situations 
the  "complexity"  of  the  definable  functions  can  be  controlled  in  a 
way  sufficient  to  ensure  continuity. 

One  example  arises  in  the  semialgebraic  category.  Fix  a  (finite) 
semialgebraic  stratification  5  of  a  compact  definable  A',  and  consider 
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definable  semialgebraic  functions  algebraic  with  respect  to  this  strati¬ 
fication  (that  is  such  that  the  restriction  of  the  function  to  any  stratum 
5  G  5  is  a  polynomial  P5).  The  resulting  linear  space  (filtered  by  the 
subspaces  of  polynomials  of  bounded  degree)  can  be  equipped  with 
the  structure  of  a  Banach  space,  by  completing  the  family  of  semi¬ 
norms  ||P||5.k  =  max5e5  ||^5||c"/  where  n  =  dinlA^  Then 
becomes  a  continuous  (non-linear)  functional  on  this  Banach  space. 

The  proof  results,  essentially,  from  the  Bezout  theorem  (mimicking 
Thom-Milnor  theory):  the  total  number  of  critical  points  graph  of 
a  polynomial  of  degree  D  on  a  fixed  semi-algebraic  set  is  bounded 
by  0{D'^).  The  generalization  to  increasing  (refined)  stratifications  is 
straightforward. 

Integration  itself  defines  a  natural  topology  for  Def(A")  on  which 
integration  is  continuous.  Define  the  6-neighborhood  of G  Def(A") 
as  the  intersection  of  the  6-neighborhood  (definable  functions  with 

6-close  graphs)  with  those  functions  g  G  Def (A")  satisfying  |/^.  /  —  g  [dyj  |  < 
6.  This  provides  a  basis  for  an  O  topology  on  Def(A").  As  a  conse¬ 
quence  of  Lemma  6,  the  definition  is  independent  of  the  use  of  [r/\ J 
or  \(l\]. 

The  interested  reader  may  speculate  on  other  function  space  topolo¬ 
gies  on  Def{A"). 

1.4.4.  Duality  and  links.  There  is  an  integral  transform  on  CF{X) 
that  is  the  analogue  of  Poincare- Verdier  duality  [37].  It  extends  seam¬ 
lessly  to  integrals  on  Def(A')  by  means  of  the  following  definition. 

Definition  17.  The  duality  operator  is  the  integral  trans¬ 
form  V  :  CF{X)  CF{X)  given  by 

(1.21)  T)h[x)  =  lim  [  IiIb  (v)d\. 

<^-0  •  ./  V 

where  is  an  open  metric  ball  of  radius  f . 

We  extend  the  definition  to  V  :  Def(A')  — >  Def(A')  by  integrating 
with  respect  to  [f/yj  or  [dy],  interchangeably,  via: 

Lemma  18.  Vh  is  well-defined  on  Def(A')  and  independent  ofivhether 
the  integration  in  (1.21)  is  with  respect  to  [d\J  or  [d.v]. 

For  a  continuous  definable  function  h  on  a  manifold  .\I,  Vh  = 

(  — 1)'*""^^/?,  as  one  can  verify  by  combining  Eqns.  [1.9]  and  [1.21]. 

This  is  commensurate  with  the  result  of  Schapira  [36]  that  V  is  an 
involution  on  CF{X). 

Theorem  19.  Duality  is  involutive  on  Def(A");  V  o  Vh  =  h. 
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One  can  define  related  integral  transforms.  For  example,  the 
LINK  of  he  CF{X)  is  defined  as 


(1.22) 


The  link  of  a  continuous  function  on  an  n-manifold  M  is  multipli¬ 
cation  by  1  +  (—1)”,  as  a  simple  computation  shows.  In  general, 
A  Id  -  P,  where  Id  is  the  identity  operator. 

1.5.  Integral  transforms.  Integration  with  respect  to  Euler  char¬ 
acteristic  over  CF{X)  has  a  well-defined  and  well-studied  class  of 
integral  transforms,  expressed  beautifully  in  Schapira's  work  on  in¬ 
version  formulae  for  the  generalized  Radon  transform  in  f/\  [35]. 
Integral  transforms  with  respect  to  [dx\  and  \(lx]  are  similarly  ap¬ 
pealing,  with  applications  to  signal  processing  as  a  primary  moti¬ 
vation.  Examples  of  interesting  definable  kernels  for  integral  trans¬ 
forms  over  Euclidean  include  ||x  — ?y||,  (x,  y),  and  g{x  —  y)  for  some 
g.  These  evoke  Bessel  (Hankel)  transforms,  Fourier  transforms,  and 
convolution  with  g  respectively.  The  choice  between  [dx\  and  \d\] 
makes  a  difference,  of  course,  but  can  be  amalgamated.  Example: 
for  fixed  kernel  K,  one  can  consider  the  mixed  integral  transform 
h  ^  /i/i  [c/xj  —  In  the  case  of  A^(x,^)  =  (x.<J),  this 

transform  takes  1,^  for  A  compact  and  convex  to  the  'width'  of  A 
projected  to  the  ^-axis. 

1.5.1.  Convolution,  On  a  vector  space  V"  (or  Lie  group,  more  gen¬ 
erally),  a  convolution  operator  with  respect  to  Euler  characteristic  is 
straightforward.  Given  f,g  e  CF{V),  one  defines 


(1.23) 


Convolution  behaves  as  expected  in  CF{V).  By  reversing  the  or¬ 
der  of  integration,  one  has  immediately  that  Jyf^gd\  —  [y  f  d\  fygd\. 
There  is  a  close  relationship  between  convolution  and  the  Minkowski 
sum,  as  observed  in,  e.g.,  [21]:  for  A  and  B  convex  and  closed  1  4  * 

1b  =  1>^+B/  cf.  [38,  36].  Convolution  is  a  commutative,  associative 
operator  providing  CF{V)  with  the  structure  of  an  (interesting  [8]) 
algebra. 

Convolution  is  well-defined  on  Def(V^)  by  integrating  with  re¬ 
spect  to  [r/xj  or  \dx].  However,  the  product  formula  for  f  f  fails  in 
general,  since  one  relies  on  the  Fubini  theorem  to  prove  it  in  CF{\'). 

1.5.2.  Linearity.  The  nonlinearily  of  the  integration  operator  pre¬ 
vents  most  straightforward  applications  of  inversion  formulae  a  la 
Schapira.  Fix  a  kernel  K  e  Def(A''  x  V")  and  consider  the  integral 
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transform  7/(:  :  Def(X)  — +  Def(y')  of  the  form  (7/^-/^)(?y)  =  li{x)K{.i\y)[dx\{.v). 

In  general,  this  operator  is  non-linear,  via  Lemma  12.  However,  some 
vestige  of  (positive)  linearity  survives  within  CF. 

Lemma  20.  The  integral  transform  Ti^-  is  positive-linear  over  C F'*'  {X)  = 
Def(X,N). 

This  implies  in  particular  that  when  one  convolves  a  function  li  6 
CF"^(R”)  with  a  smoothing  kernel  (e.g.,  a  Gaussian)  as  a  means  of 
filtering  noise  or  taking  an  average  of  neighboring  data  points,  that 
convolution  may  be  analyzed  one  step  at  a  time  (decomposing  h). 

Integral  transforms  are  not  linear  over  all  of  CF(  A"),  since  f  — /)[(/\'J  7^ 

—  f  However,  integral  transforms  which  combine  [(/\J  and 

[dy]  compensate  for  this  behavior.  Define  the  measure  [dy]  to  be  the 
average  of  [dyj  and  [dy]; 

(1.24)  ^  h[dx]  =  ^  h[dx\  +  j  dfdy]^  . 

Theorem  21.  Any  integral  transform  of  the  form 

(1.25)  {rKh){y)  =  h{x)K{.r,  y)[d,x](.r) 

is  a  linear  operator  C F {X)  — >  Def(>'). 

As  a  simple  example,  consider  the  transform  with  kernel  K{xX)  — 

(.T,  ^).  The  transform  of  l,i  with  respect  to  [dy]  for  A  compact  and 
convex  equals  a  'centroid'  of  A  along  the  c^-axis:  the  average  of  the 
maximal  and  minimal  values  of  ^  on  OA.  Note  how  the  dependence 
on  critical  values  of  the  integrand  on  DA  reflects  the  Morse-theoretic 
interpretation  of  the  integral  in  this  case. 

Integration  with  respect  to  [dy]  seems  suitable  only  for  integral 
transforms  over  CF.  On  a  continuous  integrand,  the  integral  with 
respect  to  [dy]  either  returns  zero  (cf.  the  integral  of  Rota  [34])  or 
else  the  integral  with  respect  to  [dyj,  depending  on  the  parity  of  the 
dim  A'^,  via  Eqn.  [1.18]. 

1.6.  Applications  of  definable  Euler  integration.  The  Euler  cal¬ 
culus  on  CF  is  quite  useful;  the  extension  to  Def  deepens  this  utility 
and  opens  new  potential  applications,  of  which  we  highlight  a  few. 

1.6.1.  Sensor  networks.  An  application  of  Euler  integration  over 
CF{X)  to  sensor  networks  problems  was  initiated  in  [1].  Consider  a 
space  X  whose  points  represent  target-counting  sensors  that  scan  a 
workspace  IE.  Target  detection  is  encoded  in  a  SENSING  RELATION 
S  C  IE  X  A'  where  {w,  x)  €  >S  iff  a  target  at  w  is  detected  by  a  sensor 
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at  X.  Assume  that  sensors  count  the  number  of  sensed  targets,  but 
do  not  locate  or  identify  the  targets.  The  sensor  network  therefore 
induces  a  TARGET  COUNTING  FUNCTION  h  :  A'  N  of  the  form  h  = 
Ha  where  Ua  is  the  TARGET  SUPPORT  —  the  set  of  sensors  which 
detect  target  o.  Euler  integration  allows  for  simple  enumeration; 


Theorem  22  ([!]).  Assume  h  e  CF{X)  and  \  {U„)  =  N  ^  0  for  oil 
o.  Then  the  number  of  targets  in  W  is  precisely  h(i\. 


Since  the  target  count  is  presented  as  an  integral,  it  is  possible  to 
accurately  estimate  the  answer  when  the  integrand  h  is  known  not 
on  all  of  A'  (a  continuum  of  sensors  being  an  idealization)  but  rather 
on  a  sufficiently  dense  grid  of  sample  points  (physical  sensors  in  a 
network). 

The  R-valued  theory  aids  in  establishing  expected  values  of  tar¬ 
get  counts  in  the  presence  of  confidence  measures  on  sensor  read¬ 
ings.  Let  J\f  =  {xj}  denote  a  discrete  set  of  sensor  nodes  in  R",  and 
assume  each  sensor  returns  a  target  count  h{x^)  C  N  and  a  fluctuation 
measure  c(xj)  €  [0, 1]  obtained,  say,  by  stability  of  the  reading  over  a 
time  average.  View  h  as  a  sampling  over  Af  of  the  true  target  count 
/  =  1[;„.  Assume  that  nodes  with  fluctuation  reading  0  have  per¬ 

fect  information  {h  =  /  at  x,)  and  that  c  correlates  with  error  |/  -  /)|. 
Assume  that  sensor  nodes  J\f  are  the  vertex  set  of  a  triangulation  T. 

The  integral  of  an  extension  of  /  over  a  triangulation  gives  a  ter¬ 
rible  approximation  to  f  h  d^:  an  error  of  ±1  on  K  nodes  can  cause 
a  change  in  the  integral  of  order  K.  More  specifically,  if  /;  —  /  }  v, 
where  e  :  A/”  — »  {  — 1, 0, 1}  is  an  error  function  that  is  nonzero  on  a 
sparse  subset  M'  C  M ,  then,  for  certain  infelicitous  choices  of  Af , 

\!h-ss\  =  m- 

A  M-valued  relaxation  can  mitigate  errors  by  using  fluctuation 
c  as  a  weight.  Let  N{i)  be  a  collection  of  neighboring  nodes  to 
where  neighborhood  can  be  defined  via  distance  (if  available)  or 
edge*distance  (in  an  ad  hoc  network  or  triangulation).  Define  //  to 
be  the  result  of  averaging  the  value  at  j;,  C  M  over  with  r  as  a 
weight.  Specifically, 


(1.26) 


r,  ,  JlyeN{,)iiyMy) 


This  nearest-neighbor  convolution  damps  out  local  variations.  The 
resulting  integral  with  respect  to  [dx\  will  tend  to  mitigate  localized 
errors,  thanks  to  the  Morse-theoretic  formula. More  numerical  inves¬ 
tigation  is  warranted. 
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Such  averaging  leads  to  non  integer-valued  integrands.  By  us¬ 
ing  integration  with  respect  to  [c/yj  or  [d\]  for  upper/lower  semi- 
continuous  integrands  associated  to  an  averaged  signal  h,  one  ob¬ 
tains  an  expected  value  of  /  h  dx-  This  can  be  particularly  illuminat¬ 
ing  when  a  network  has  incomplete  information,  e.g.,  a  hole.  Holes 
in  a  network  can  be  modeled  by  setting  the  confidence  measure  c  to 
zero  and  averaging. 

1.6.2.  Numerical  integration.  Though  integration  with  respect  to 
Euler  characteristic  has  a  lengthy  history,  there  appears  to  be  no 
treatment  of  numerical  integration,  even  in  the  simpler  setting  of 
C'F'^(R”).  The  central  problem  (in  the  constmctible  and  definable 
categories)  is  how  to  estimate  h  given  the  values  of  h  on  a  discrete 
subset  of  X.  As  in  the  case  of  numerical  integration  for  Riemann 
integrals,  one  typically  assumes  something  about  the  features  of  li 
and/ or  the  density  and  extent  of  the  sampling  set.  In  [1],  the  present 
authors  give  a  formula  for  estimating  J  lid\  given  a  discrete  sam¬ 
pling  of  /i  €  C'F(R^)  which  correctly  samples  connectivity  data  of 
excursion  sets.  This  formula  generalizes  to  the  definable  category: 

Proposition  23.  For  h  G  Def(R^)  continuous,  J  h  [dxj  = 


(1.27)  /  pQ{h  >  .s}  +  /^{/^.  >  — .s}  —  do{/i  <  .s}  —  ds. 


where  do(')  =  dim  the  zeroth  Betti  number. 

The  value  of  Proposition  23  is  that  it  allows  for  computation  based 
on  do  quantities.  Such  connectivity  data  are  easily  obtained  from  a 
discrete  sampling  via  clustering.  We  have  implemented  this  formula 
in  software.However,  for  more  general  integration  domains  than  R“, 
duality  formulae  are  less  helpful.  One  general  result  on  refinement 
follows  from  continuity  of  the  integral  operator. 

Theorem  24.  For  h  g  Def(A')  continuous,  let  hpi,  be  the  pieceudse- 
linear  function  obtained  from  sampling  h  on  the  vertex  set  of  a  triangulation 
T  of  X.  As  the  sampling  and  triangulation  are  refined. 


(1.28) 


This  result  relies  crucially  on  continuity  and  does  not  apply  to 
CF{X).  A  more  desirable  result  would  be  a  measure  of  how  far  a 
given  sampling  is  from  the  true  integral.  This  seems  challenging. 
We  note  that  the  Morse-theoretic  formulae  [1.16]-[1.17]  allow  one  to 
reduce  the  domain  of  an  integral  to  a  (typically  finite)  set  of  critical 
points.  This  'focusing'  property  of  integration  over  Def  should  be  a 
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starting  point  for  good  numerical  algorithms,  especially  for  integral 
transforms. 


2.  Unimodal  Category 

2.1.  Introduction.  Here  we  introduce  a  novel  topological  approach 

to  the  problem  of  decomposing  distributions  /  :  — >  [0,  oc)  into  a 

convex  combination  of  basis  distributions.  Instead  of  employing  a 
decomposition  into  analytically  defined  {e.g.,  normal)  distributions, 
we  propose  a  decomposition  into  topologically  defined  factors.  Specif¬ 
ically,  we  consider  the  decomposition  of  a  distribution  into  a  sum  of 
UNIMODAL  distributions:  those  with  a  single  maximum  value  and 
no  other  extrema.  Such  a  decomposition  is  not  uniquely  defined; 
however,  the  minimal  number  of  unimodal  summands  is.  This  UNI¬ 
MODAL  CATEGORY,  UCAT(/),  is  a  coarse  measure  of  complexity  for  a 
distribution.  The  computation  of  UCAT  is  not  obvious. 

2.1.1.  Motivation.  There  are  several  contexts  within  which  the  ques¬ 
tion  about  the  number  of  summands  in  a  unimodal  decomposition 
arises.  The  first  such  context  is,  obviously,  the  statistical  one.  Indeed, 
unimodal  distributions  are  the  primal  building  blocks  of  statistical 
models.  Essentially  all  classical  probability  distributions,  including 
normal,  Poisson,  Gamma,  Beta,  Bernoulli,  and  more,  are  all  uni¬ 
modal.  The  methodology  of  statistical  modeling  essentially  forces 
one  to  assume  that  the  presence  of  several  modes  in  a  distribution  is 

a  consequence  of  its  being  a  mixture  of  several  distributions,  and  the 
relationship  between  the  number  of  modes  (essentially,  the  number 
of  local  maxima  of  the  density,  in  the  multivariate  case)  and  the  num¬ 
ber  of  ''components"  of  the  mixture,  i.e.,  the  number  of  summands 
in  the  convex  decomposition,  has  been  studied  by  many  authors. 

To  justify  the  radical  difference  of  our  setup  from  the  traditional 
statistical  one  —  our  lack  of  any  assumptions  about  the  structure  of 
the  summands  beyond  unimodality  —  we  can  invoke  two  consider¬ 
ations: 

(1)  The  very  variety  of  unimodal  building  blocks  in  standard 
statistical  models  suggests  that  one  should  try  to  abstract 
away  any  specific  distribution,  retaining  only  the  minimal 
topological  properties. 

(2)  Any  specific  analytic  form  of  a  density  binds  the  distribution 
to  some  fixed  coordinate  system.  Adopting  a  completely 
coordinate-free  decomposition  will  inexorably  relax  to  a  topo¬ 
logical  approach. 
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We  envisage  applications  of  this  topological  decomposition  to  a 
variety  of  contexts  within  statistics,  as  well  as  data  analysis  and  visu¬ 
alization  (where  decomposition  can  play  a  role  in  efficient  encoding 
of  the  images  or  of  multidimensional  data). 

The  particular  application  inspiring  this  paper  comes  from  sig¬ 
nal  analysis  and  reconstruction  in  a  problem  of  target  enumeration 
by  means  of  a  field  of  sensors  [?].  In  that  setting,  one  has  a  finite 
number  of  targets  located  throughout  a  domain  D.  The  domain  is 
assumed  to  be  filled  with  a  network  of  primitive  sensors.  A  sensor 
atx  E  D  returns  a  number  h{x)  6  N  which  represents  a  count  of  how 
many  targets  it  senses.  The  problem  is  to  compute  the  number  of  tar¬ 
gets  given  the  function  h  :  D  ^  N.  Assume  each  target  o  has  a  com¬ 
pact  contractible  support  U„  on  which  it  is  observed  by  all  sensors 
within.  Then,  the  paper  [?]  computes  the  number  of  targets  based 
on  the  Euler  characteristics  of  excursion  sets  of  h.  When,  however, 
one  assumes  a  sensor  modality  which  reads  a  [0,  oo)-valued  signal 
instead  of  an  integral  count,  the  methods  of  [?]  do  not  apply,  and  the 
techniques  of  this  report  take  over.  For  example,  in  a  network  with 
acoustic  or  infrared  sensors,  it  is  more  reasonable  to  expect  sensors 
to  return  a  real-valued  signal  which  sums  the  contributions  of  in¬ 
range  targets,  these  contributions  varying  in  intensity  as  a  function 
of  range-to-target.  Assuming  a  distribution  of  target  sensory  impact 
is  a  very  weak  assumption;  however,  the  methods  of  this  paper  yield 
lower  bounds  on  the  enumeration  problem. 

2.1.2.  Statement  of  Results.  The  extent  of  this  report  is  to  define 
the  unimodal  category,  give  a  complete  characterization  and  method 
of  computation  in  the  univariate  case,  and  provide  some  key  steps  to 
understanding  decompositions  of  higher-dimensional  distributions. 
The  specific  contributions  are  as  follow: 


(1)  We  observe  (Lemma  30)  that  the  unimodal  category  is  in¬ 
variant  under  the  right-action  of  a  homeomorphism  on  the 
domain.  This  demonstrates  the  topological  naturality  of  our 
definitions. 

(2)  We  give  a  complete  description  (Theorem  36)  of  the  uni¬ 
modal  category  in  the  univariate  case,  along  with  a  simple 
algorithm  for  its  computation. 

(3)  We  demonstrate  that  the  univariate  and  bivariate  cases,  the 
unimodal  number  is  a  function  of  the  Reeb  graph  of  the  dis¬ 
tribution,  labeled  by  critical  values  (Proposition  39).  This 
implies  that  the  unimodal  category  is  purely  a  function  of 
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Morse-theoretic  data,  providing  a  discretization  of  the  prob¬ 
lem. 

(4)  We  prove  a  result  (Proposition  41)  about  the  'persistence'  of 
unimodal  category  with  respect  to  translating  the  graph  of 
the  distribution.  Shifting  vertically  can  only  decrease  the 
unimodal  category,  and  this  number  stabilizes  to  the  value 
two. 

(5)  For  distributions  in  arbitrary  dimensions,  we  show  (Proposi¬ 
tion  ??)  that  a  minimal  unimodal  decomposition  exists  with 
"depth"  (the  degree  of  overlaps  of  supports)  bounded  above 
by  the  dimension  of  the  domain.  We  believe  this  result  will 
be  useful  in  computing  the  unimodal  category  by  means  of 
topological  tools. 

2.1.3.  Related  work.  Morse  structures  associated  to  mixtures  of 
multivariate  normal  distributions  are  discussed  in  [32].  This  ques¬ 
tion  already  seems  to  be  of  interest  to  statisticians  —  even  the  simple 
univariate  case  is  discussed  in  detail  across  several  papers  [16,  4,  33], 
while  the  mixture  of  non-Gaussian  unimodal  densities  is  considered 
in  [24,  26].  In  particular,  it  is  known  that  in  mixtures  of  normal  uni¬ 
variate  distributions,  the  number  of  modes  cannot  exceed  the  num¬ 
ber  of  components,  a  result  which  does  not  hold  in  higher  dimen¬ 
sions  (compare  [11]). 

The  questions  posed  in  this  report  are  mostly  complementary  to 
the  recent  works  [9,  10,  12,  ?,  30]  considering  the  estimation  of  the 
topology  of  data  sets  generated  from  noisy  point-clouds.  This  is  one 
set  of  works  demonstrating  the  increasing  interest  in  how  topologi¬ 
cal  and  statistical  phenomena  entwine. 

2.2.  Unimodal  category.  The  following  definitions  are  common. 
We  follow  the  notation  of  [14]. 

Definition  25.  For  A'  a  topological  space,  the  Lusternik-Schni- 
RELMANN  CATEGORY  of  X,  LSCAT(A''),  is  the  minimum  number  of 
open  sets  contractible  in  X  which  cover  A'.  The  GEOMETRIC  CATE¬ 
GORY  of  A",  gcat(A"),  is  the  minimum  number  of  open  sets  homotopic 
to  a  point  which  cover  A'. 

Some  authors  (including  those  of  [14])  use  a  reduced  category, 
which  measures  the  minimal  number  of  open  sets  minus  one.  We  do 
not  follow  this  convention.  Geometric  category  gcat  is  a  homeomor- 
phism  invariant  of  a  space,  and  L-S  category  LSCAT  is  a  homotopy 
invariant.  There  are  numerous  deep  connections  between  category 
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and  critical  point  theory  (the  classical  motivation  for  the  subject),  dy¬ 
namical  systems,  homotopy  theory,  and  symplectic  topology. 

We  introduce  a  new  variant  of  category  for  distributions  based 
on  decomposition  into  unimodal  factors. 

Definition  26.  For  A'  a  topological  space,  let  D  =  2)(A')  denote 
the  set  of  all  compactly  supported  continuous  distributions  /  :  A' 
[0,oo). 

Definition  27.  A  distribution  G  2)  is  said  to  be  UNIMODAL  if 
the  upper  excursion  sets  oo))  have  the  homotopy  type  of 

a  point  for  all  0  <  c  <  M  and  are  empty  for  all  r  >  M.  Such  a  ii  has 
M  as  its  maximal  value. 

We  will  refer  to  the  nonempty  upper  excursion  sets  C  A  as 
being  CONTRACTIBLE,  though  it  must  be  clarified  that  such  sets  are 
contractible  in  themselves  as  opposed  to  contractible  in  A'.  The  lat¬ 
ter  would  be  more  in  line  with  the  definitions  used  in  Lustemik- 
Schnirelmann  theory,  but  would  render  the  theory  useless  for  most 
applications  (where  A'  =  R"). 

Definition  28.  Fix  a  norm  =  H-H  on  R^**.  The  UNIMODAL  I'- 
CATEGORY  of  a  distribution  /  C  2) (A")  is  defined  as  the  minimal 

number  UCAT"  of  unimodal  distributions  =  1 . UCAT''  on  X 

such  that  /  is  pointwise  the  z^-norm  of  the  collection  («o)-  Specifi¬ 
cally,  f{x)  =  ||(zio(x))||  for  all  x  C  A'. 

The  most  natural  and  fundamental  example  is  the  minimal  num¬ 
ber  of  unimodal  distributions  required  to  represent  a  given  distribu¬ 
tion  as  a  sum  of  unimodals.  Summation  of  the  components  corre¬ 
sponds  to  the  1-norm  on  vectors,  leading  to  the  following  general¬ 
ization. 

Example  29.  The  unimodal  ^-category  of  a  distribution  /  c 

2)  is  the  minimal  number  of  unimodal  distributions  ?/,).  o  =  1 . UCAT'' 

such  that  /  is  pointwise  an  combination.  Specifically, 


(2.1) 


p  =  oo  :  f{x)  =  max  {ztft(a:)} 


(2.2) 


Lemma  30.  Any  unimodal  u-category  is  invariant  under  the  right- 
action  of  the  homeomorphisrn  group. 

2.3.  Unimodal  1-category.  Motivations,  etc. 
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General  observations. 

Remark  31.  A  trivial  upper  bound  on  UCAT*(/)  is  the  number  of 
local  maxima  of  /.  The  corresponding  unimodal  summands  can  be 
found  supported  on  a  small  neighborhood  of  the  basins  of  attraction 
of  the  gradient  flow. 

Definition  32.  Assume  that  a  unimodal  decomposition  of  /  = 
is  given.  A  set  U  C  sup(/)  is  called  MAX-FREE  if  U  does  not 
contain  any  of  the  critical  points  of  any  /„.  For  any  max-free  U  we 
denote  by  h{U)  the  DEPTH  of  U,  the  number  of  functions  of  the  uni¬ 
modal  decomposition  not  vanishing  on  U: 

h{U)  =  #{a  :  f/nsiip(a«)  0). 

The  following  trivial  lemma  is  useful: 

Lemma  33.  The  depth  of  an  open  max-free  set  U  is  bounded  from  below 

by 

max^  / 
rnaxac;  / 

Solving  the  problem  of  computing  unimodal  1-category  is  useful 
in  other  contexts; 

Lemma  34.  Tor  all  f  e  2),  UCAT''(/)  =  UCAT*  (/*/''). 

Univariate  distributions.  We  commence  with  a  computation  of 
the  unimodal  1-category  on  2)(R').  We  assume  that  this  function 
has  isolated  critical  points.  Up  to  a  homeomorphism,  this  distribu¬ 
tion  is  completely  characterized  by  the  up-down  sequence  of  critical 
values  corresponding  to  local  minima  and  maxima, 

0  =  no  <  n?o  >  7?!  <  . . .  <  nik  >  n*.  =  0; 

where  n,  =  /(.C2,);  ni,  =  f{x2i+\);  .rn  <  Xi  <  ...  <  x-ir;  and  i  = 
0, ....  A,  counting  the  initial  and  final  point  of  the  support  as  local 
minima. 

Proposition  35.  If  an  open  interval  (j‘2,,  xoj)  bounded  by  local  min¬ 
ima  is  max-free  (for  some  unimodal  decomposition  of  f),  then 

(2.3)  n,  -  m,  +  n,+  i  -  . . .  -  -F  n^  >  0. 

Consider  the  open  intervals  with  endpoints  at  the  local  minima 
{i2i}i-  Call  such  an  interval  FORCED-MAX  if  the  inequality  (2.3)  is 
violated  there.  Obviously,  forced  max  intervals  form  an  ideal:  any 
interval  containing  a  forced-max  interval  is  itself  forced-max. 
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Theorem  36.  Let  f  G  2)(R)  have  maximal  values  {m^)^  and  minimal 
values  {ui)^  ordered  according  to  the  critical  point  order  in  the  domain. 

Then  UCAT*(/)  is  equal  to  the  maximal  number  of  non-intersecting  forced- 
max  intervals: 

(2.4) 

UCAT^(/)  =  inax{A'  :  {x2^o,X2i,),  {x2^,.X2^2), forced  max}  . 

It  is  clear  that  for  any  collection  of  N  nonintersecting  forced-max 
intervals,  the  number  of  summands  in  the  unimodal  decomposition 
cannot  be  less  than  N :  each  forced-max  interval,  trivially,  contains 
a  critical  point  of  at  least  one  of  the  functions  of  the  decomposition. 

The  following  algorithm  yields  an  explicit  unimodal  decomposition 
and,  simultaneously,  a  collection  of  nonintersecting  forced-max  in¬ 
tervals,  one  for  each  summand. 

We  construct  the  functions  Ua  iteratively,  left  to  right,  according 
to  the  algorithm  Sweep.  This  entails  sweeping  /  from  the  left  and 
pulling  out  unimodal  factors  which,  on  their  descent,  compensate 
for  the  remaining  factors  as  much  as  possible  by  descending  accord¬ 
ing  to  the  (positive)  slope  df.  Here,  the  functions  are  the  differ¬ 
ences  of  heights  of  the  interleaving  curves,  which  have  the  same  dif¬ 
ferential  as  /  on  the  intervals  where  /  is  nonincreasing.  A  new  sum¬ 
mand  starts  after  the  previous  curve  hits  the  x  axis. 


Algorithm  1  {ua]  =  Sweep(/) 

Require;  /  E  'D{R)  with  minima  rii  —  f{j  •2,)S  and  maxima  //i,  ~ 

f{X2i-l)l 

1;  ifo  0  ;  o  4=  1  ;  4=  / 

2;  while  g„  7^  0  do 

3:  y„  4=  first  maximum  of  from  left 

4:  Ua\  {-OO^lJa]  Oa 

5:  dua\  iVe,,  00)  4=  inm{df,  0) 

6:  Ua  inax(7i^j,  0) 

7:  increment  a 

9:  end  while 
10:  return 


The  proof  that  this  construction  is  minimal  follows  from  the  ob¬ 
servation  that  the  local  minima  spanning  the  curves  in  the  graphical 
construction  form  a  max-forced  partition  of  the  support  of  /. 

We  have  observed  that,  in  the  univariate  case,  the  unimodal  1- 
category  is  a  function  of  the  critical  values  and  the  order  in  which 
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they  appear;  however,  there  is  not  a  strict  dependence  on  this  or¬ 
dering.  For  example,  in  accordance  with  Lemma  30,  UCAT^  must  be 
invariant  under  reversing  the  order  of  the  critical  points  (.r  ►  — x). 

There  is  in  fact  a  weaker  dependence  on  the  labeled  Reeb  graph  of 
the  distribution. 

Definition  37.  The  Reeb  graph  of  a  Morse  function  /  on  a 
smooth  manifold  M  is  the  quotient  space  of  M  with  respect  to  the 
relation  declaring  two  points  equivalent  if  they  belong  to  the  same 
connected  component  of  a  level  set  of  /.  In  particular,  the  vertices  of 
the  graph  are  given  by  the  connected  singular  level  sets  of  /. 

Remark  38.  In  the  univariate  case,  the  Reeb  graph  is  a  rooted 
metric  tree,  i.e.,  a  rooted  tree  with  lengths  attached  to  its  edges.  The 
construction  of  this  tree  from  an  excursion  is  well-known  in  combi¬ 
natorics  and  probability  theory  (cf.  [31]).  From  our  construction  it 
follows  that  for  a  univariate  distribution  UCAT^  depends  only  on  its 
Reeb  graph.  In  other  words,  one  can  swap  the  subexcursions  bor¬ 
dering  at  a  local  maxima  without  affecting  the  unimodal  number. 

2.3.1.  Toward  uiiiHivariate  decompositions.  Multivariate  distributions 
introduce  a  number  of  complexities.  We  sketch  a  few  results  in  this 
section,  reserving  a  more  complete  treatment  for  a  later  work. 

An  indication  of  the  troubles  arising  in  the  multivariate  case,  we 
notice  that  an  analogue  of  Proposition  35  is  not  valid  there. 

The  dependence  of  the  unimodal  category  on  the  Reeb  graph 
noted  in  Remark  38  can  be  extended  to  planar  distributions: 

Proposition  39.  The  imimoda}  l-category  of  f  e  T){W)for  n  — 
is  a  function  of  the  combinatorial  type  of  the  Reeb  graph  of  f  labeled  by 
critical  values. 

Remark  40.  In  higher  dimensions,  the  Reeb  graph  does  not  carry 
enough  information:  for  example,  a  critical  point  of  signature  (2, 1) 
is  a  vertex  of  valence  2  on  the  Reeb  graph,  with  no  information  about 
the  change  of  topology  under  the  surgery  defined  by  the  point.  In  di¬ 
mension  2,  on  the  contrary,  the  change  of  topology  is  unambiguous. 

Unimodal  1 -category  is  not  solely  a  function  of  the  critical  point 
and  Morse  index  data  of  /,  but  rather  depends  upon  the  critical  val¬ 
ues  as  well.  This  dependence  can  be  viewed  in  the  following  param¬ 
eterized  sense. 

Proposition  41.  Fix  a  (non-nnimodal)  C*  distribution  f  e  D(M")- 
For  any  0  <  p  <  oo,  the  unimodal  p-category  of  the  shifted  distribution 
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obtained  by  adding  Cls,ip/  to  f  for  constant  C  >  0  is  a  non-increasing 
function  of  C  which  stabilizes  to 

(2.5)  lim  UCAT''  (/  +  Clsu,./)  =  1  +  UCAT"(/)  =  1  +  gcat(sup  /). 

C— »oo 

If  one  normalizes  the  distributions  /  +  Clgup/  to  have  unit  mass, 
then  one  sees  clearly  the  effect  of  increasing  C  is  to  reduce  the  total 
variation. 

2.4.  Unimodal  oo-category.  Motivations,  etc. 

General  observations.  Computing  the  unimodal  oo-category  in 
the  univariate  setting  is  trivial. 

Lemma  42.  For  any  f  e  S)(IR*),  UCAT‘^(/)  equals  the  number  of  local 
maxima  of  f. 

Lemma  43.  UCAT'^  is  itwariant  under  the  left  action  o/Homeo[0.  oc). 

2.5.  Bounds  on  unimodal  category.  In  cases  where  exact  com¬ 
putation  of  unimodal  category  is  difficult,  one  turns  to  bounds  based 
on  analytic,  geometric,  or  topological  features  of  the  distribution. 

2.5.1.  Unimodal  0-category.  The  monotonicity  result  allows  one  to 
define  a  unimodal  category  for 

Definition  44.  The  unimodal  0-category,  UCAI^,  is  defined  as  the 
limit  of  UCATP  as  p  ^  0. 

This  limit  is  well-defined. 

Theorem  45.  For  any  /  e  S,  UCAT“(/)  =  gcat(sup(/)),  the  geo¬ 
metric  Lusternik-Schnirelmann  category  of  the  support  of  f. 

Proof.  From  Lemma  34,  UCAT'(/)  =  UCAT*  As  p  ^  0,  f'>‘ 

approaches  l.supi/)- 
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